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Abstract 

The  Mercury  generator  is  a  ~6  MV,  -360  kA,  -50  ns 
magnetically-insulated  voltage  adder  that  recently  became 
operational  at  the  Naval  Research  Laboratory  [1,2], 
providing  a  radiation  environment  of  interest  to  both  the 
DOE  and  DoD  communities.  To  establish  the 
bremsstrahlung  dose  and  dose-rate  production  of 
Mercury,  a  series  of  large-area  electron-beam  (e-beam) 
diode  experiments  are  performed  in  which  the  e-beam  is 
deposited  onto  either  a  carbon  or  an  optimized-thickness 
tantalum  converter  package.  In  this  paper,  experimental 
results  are  presented  and  compared  with  predictions  from 
LSP  (PIC)  simulations  [3]  and  ITS/Cyltran  (Monte  Carlo) 
calculations  [4]. 


I.  INTRODUCTION 
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Figure  1.  Schematic  of  the  front-end  of  the  Mercury  pulsed 
power  generator,  showing  current  monitors,  blade  e-beam 
diode,  vacuum  vessel  materials,  and  radiation  diagnostics. 


To  quantify  the  radiation  environment  of  the  Mercury 
generator,  experiments  are  performed  utilizing  both 
carbon  (C)  and  tantalum  (Ta)  converters.  An  array  of 
CaF2  thermoluminescent  dosimeters  (TLDs)  and 
scintillator  photodiodes  (PDs)  are  used  to  measure  the 
dose  and  dose-rate  at  various  angles  and  distances  from 
the  source.  For  a  particular  bremsstrahlung  converter 
package  (including  the  anode  foil,  converter  material,  and 
electron  beam  stop),  the  voltage-dependent  electron 
dynamics  in  the  anode-cathode  gap  directly  affect  the 
magnitude  and  angular  dependence  of  the  radiation.  A 
series  of  steady-state  PIC  simulations  are  performed  using 
LSP  [3]  to  model  the  electron  distribution  incident  on  the 
anode  throughout  the  Mercury  voltage  pulse.  These 
simulations  include  both  power-flow  electrons  [5,6]  born 
in  the  magnetically-insulated  inductive  voltage  adder 
(MIVA)  and  electrons  born  in  the  diode.  Electron 
scattering  and  energy  loss  in  the  bremsstrahlung  converter 
package  are  modeled  self-consistently  by  LSP. 
Predictions  of  the  bremsstrahlung  production  (angular- 
dependent  dose  and  dose-rate)  are  generated  by  using  the 


PIC  electron  distributions  as  initial  conditions  for  Monte 
Carlo  radiation-transport  simulations  with  the  ITS-Cyltran 
code  [4].  These  predictions  are  compared  with 
experiments  on  Mercury  for  both  C  and  Ta  converters, 
establishing  the  bremsstrahlung  capabilities  of  the 
generator  and  benchmarking  our  predictive  capability. 


II.  EXPERIMENTAL  SET-UP 

A  schematic  of  the  front-end  of  Mercury  is  shown  in 
Fig.  1.  The  bremsstrahlung  experiments  described  here 
use  a  cylindrical  electron-beam  diode  in  which  high- 
energy  electrons  from  a  large-diameter,  annular  cathode 
impinge  upon  a  planar  anode.  To  minimize  e-beam 
backscattering  and  subsequent  ion  formation,  a  low 
atomic-material  was  used  as  the  anode  ground  plane.  On 
Ta  shots,  an  aluminum  (Al)  was  used;  for  C  shots,  the 
carbon  itself  served  this  purpose.  These  designs  inhibit 
pinching  and  result  in  an  annular  radiation  source. 
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Figure  2.  Anode  and  cathode  currents  measured  near  the  diode  on  Shots  302  (carbon  anode)  and  317  (tantalum  anode) 
are  plotted  versus  time.  The  voltage  is  calculated  using  the  Mendel  formula,  and  peaks  around  6  MV.  Forward-directed 
dose-rates  dD/dt  are  shown  in  units  of  [rad(CaF2)/ns]. 


For  Shots  280  -  306,  a  cylindrical-diode  with  an  11.55 
cm  radius  cathode  deposits  an  e-beam  onto  a  1.27-cm 
thick  C  anode,  bolted  onto  a  2-cm  thick  A1  beam  stop. 
Throughout  this  shot  series,  the  Marx  charge  and  diode 
gap  is  varied  from  50  to  75  kV  and  from  4  to  12  cm, 
respectively,  to  calibrate  diagnostics  and  to  study  MITL 
powerflow  and  diode  behavior  [2,7].  A  detailed  analysis 
of  Shot  302  is  presented  in  this  paper.  The  diode  voltage 
on  Shot  302  exceeds  6  MV  and  the  MITL  runs  “line- 
limited”  [7]  with  respect  to  the  12-cm  A-K  gap  in  the  e- 
beam  diode.  The  radiation  production  is  also 
characterized  with  a  full  array  of  x-ray  diagnostics. 

For  Shots  307  -  317,  a  cylindrical-diode  with  an  11.55 
cm  radius  cathode  deposits  an  e-beam  onto  a  Ta  converter 
package  across  a  12  cm  gap.  The  Ta  package  consists  of 
a  thin  anode  A1  foil  (with  aerodag),  multiple  foils  of  25- 
|im  Ta  as  the  bremsstrahlung  converter,  and  an  A1  beam 
stop.  The  material  thicknesses  are  varied  during  this 
series  to  provide  data  for  benchmarking  to  code 
predictions.  Experiments  and  simulations  at  6  MV 
suggest  the  on-axis  dose  is  optimized  for  75-  to  1 00-um 
of  Ta.  A  detailed  analysis  of  Shot  317  is  presented  in  this 
paper.  For  this  shot,  the  diode  voltage  is  nearly  6  MV,  the 
MITL  flow  is  “line-limited”  with  respect  to  the  12-cm 
A-K  gap,  and  the  radiation  production  is  characterized 
with  a  full  array  of  x-ray  diagnostics.  For  Shot  317,  the 
Ta  package  consists  of  a  1 .75 -pm  thick  A1  anode  foil 
(with  aerodag)  and  four  25 -[tm  thick  sheets  of  Ta, 
separated  from  a  2-cm  thick  A1  beam  stop  by  a  0.635-cm 
vacuum  gap. 

For  Shots  302  (C)  and  317  (Ta),  the  total  thickness  of 
the  bremsstrahlung  converter  package  is  about  twice  the 
CSDA  electron  range  at  6  MeV,  primarily  due  to  the  2-cm 
thick  A1  beam  stop.  The  measured  doses  could  be 
increased  by  reducing  this  beam-stop  thickness. 

To  measure  the  dose  and  dose-rate  during  the 
experiments,  an  array  of  TLDs  and  PDs  are  located  3  m 
from  the  anode  origin  at  0°,  ±20°,  and  50°  (see  Fig.  1). 
Each  PD  consists  of  a  2.5-cm  thick  plastic  scintillator  (2- 
ns  time  response)  coupled  to  a  photodiode  (0.5-ns  rise 


time)  and  is  mounted  in  a  lead  (Pb)  shield  to  view  the 
radiation  source  through  a  1.3-cm  diameter  aperture.  For 
Shot  317,  a  1.905-cm  thick  Pb  filter  is  used  in  front  of  the 
PDs  to  prevent  saturation.  Two  TLDs  are  mounted 
approximately  2.54-cm  from  each  of  the  PDs  to  measure 
the  dose  locally.  The  1x1x6  mm  CaF2  TLDs  are 
equilibrated  with  an  A1  cylinder  of  4-mm  wall  thickness; 
this  thickness  is  sufficient  to  isolate  the  TLDs  from 
externally  produced  Compton  electrons.  The  measured 
PD  signals  are  scaled  so  that  their  time-integrals  equal  the 
dose  measured  by  the  co-located  TLDs.  Therefore,  the 
PD-TLD  configuration  gives  a  time-resolved  measure  of 
dose-rate  during  the  50  ns  pulse.  Calculations  shown 
below  demonstrate  that  this  approach  is  reasonably 
accurate. 


II.  EXPERIMENTAL  RESULTS 


Anode  (Ia)  and  cathode  (Ik)  diode  currents  measured  in 
Shots  302  and  317  are  shown  in  Fig.  2.  The  current  probe 
locations  are  indicated  in  Fig.  1.  Because  the  diode 
impedance  is  large  compared  to  the  line-limited 
impedance  of  the  MITL  (ZdodB  ~  33  Q  >  ZMTL  =  23  Q),  the 
Mendel  voltage  formula  [8] 


V  =  Z 


2 _ i2  _  me 

a  k  2e 


-1 


(1) 


v "  / 

is  an  accurate  estimate  of  the  diode  voltage.  In  Eq.  (1),  the 
vacuum  impedance  Zo  =  30  Q  and  mc2/e  =  511  kV.  The 
calculated  diode  voltage  VMendd  peaks  at  approximately  6 
MV  for  both  shots.  The  measured  anode  currents  and 
voltages  are  used  in  Eq.  (3)  below  to  calculate  doses  and 
dose-rates  for  comparison  with  measurements. 

The  on-axis  dose  at  1  m  through  3.6  cm  of  A1 
(including  the  beam-stop)  is  about  87  rad(CaF,)  for  Shot 
302  and  445  rad(CaF,)  for  Shot  317.  The  peak  dose-rates 
are  2.1x1 0'  rad(CaF,)/sec  in  the  40-ns  pulse  for  Shot  302 
and  1.2xl010  rad(CaF,)/sec  in  the  36-ns  pulse  for  Shot 
317.  The  measured  on-axis  dose-rates  dD/dt,  shown  in 
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Fig.  2  with  units  of  10“  rad(CaF,)/sec,  are  accurately  fit  by 
180  IV2  for  Shot  302  and  690  IV25  for  Shot  317,  with 
currents  in  [A]  and  voltages  in  [MV]. 


III.  MODELING 


A  series  of  steady-state  LSP  [3,5,6]  simulations  are 
performed  at  several  voltages  Vk  to  assess  the  radiation 
production  from  a  given  converter  package,  i.e.  dose- 
efficiency  D(Vk,0)/Q  as  a  function  of  Vk.  Dose-efficiency 
is  defined  here  as  the  dose  D  (energy  absorbed  per  gram) 
per  deposited  electron  charge  Q  in  the  bremsstrahlung 
diode.  A  symbolic  representation  of  the  algorithm  used  in 
this  analysis  is  shown  in  Eq.  (2). 

Vk  =>  fe(x,p;Vk)  =>  y(0,E,Q; Vk )  =>  %;Vk)  (2) 

LSP  ITS  |l(E)U 


PIC  electrons  incident  on  the  anode  are  collected  in 
each  of  these  steady-state  LSP  simulations  and  are  used  to 
approximate  the  true  electron  distribution  function  (f )  at  a 
point  in  time  during  a  6-MV,  50-ns  FWHM  trapezoidal 
voltage  pulse  on  Mercury.  Each  electron  population 
collected  is  then  passed  into  a  2d  cylindrical  Monte-Carlo 
electron/photon  transport  calculation  using  Cyltran  [4], 
This  step  accounts  for  internal  absorption  of 
bremsstrahlung  photons  and  generates  an  angularly- 
dependent  point-source  spectrum  y(9,E,Q)  per  incident 
electron  charge  Q.  This  spectrum  is  then  used  in 
conjunction  with  tabulated  attenuation  and  absorption 
coefficients  [NIST]  to  generate  the  dose-efficiency 
relations  D/Q  for  each  voltage,  Vk,  and  angle,  0.  These 
D/Q  values  from  a  series  of  steady-state  voltage 
simulations  ranging  from  1  to  7  MV  are  fit  with  a  power- 
law  of  the  form  aVn  and  listed  in  Table  1. 

To  obtain  an  estimate  of  the  dose,  the  appropriate  D/Q 
power-law  is  integrated  in  time  with  the  measured  anode 
current 


D(0)  =  Jla(t) 


§[V(.),0] 


dt 


(3) 


The  estimated  dose-rate  is  simply  the  integrand  of  Eq.  (3), 
dD/dt  =  Ia(t)  (D/Q).  The  power-laws  expressed  in  Table  1 
are  appropriate  only  for  the  experimental  conditions 
modeled  here;  any  change  in  MITL  powerflow,  diode 
configuration,  converter  material  thicknesses,  TLD 
equilibration,  and  photo-diode  filtration  will  affect  the 
results.  However,  these  relations  do  suggest  qualitatively 
how  one  might  increase  the  dose  by  increasing  the  voltage 
and  serve  as  a  baseline  for  estimating  the  dose  at  any 
radial  distance  and  angle  from  a  point-source  on  this 
machine  in  this  configuration. 

Analysis  of  the  PIC  electrons  and  current  contours 
indicate  that  the  e-beam  was  incident  on  the  converter 
package  at  an  average  angle  of  20°  to  30°  near  peak 
voltage.  Qualitatively,  the  x-ray  photographs  (not  shown 
here)  indicate  that  the  e-beam  was  azimuthally-uniform 
and  justify  the  use  of  a  cylindrical  2d  model. 


A  comparison  of  measured  and  calculated  doses  and 
dose-rates  is  shown  in  Fig.  3.  In  the  0°  and  ±20° 
directions,  simulations  over-predict  the  dose  by  about 
15%.  On  the  other  hand,  predictions  of  dose-rates  at  50° 
agree  well  with  measurements.  These  comparisons 
suggest  that  the  electron  angles-of-incidence  may  not  be 
modeled  well.  For  Shot  317,  the  predicted  curves  agree 
with  measurements  during  the  rising  portion  of  the 
radiation  pulse  but  diverge  later,  suggesting  that  the 
voltage  scaling  changes  in  time. 

Table  1.  A  comparison  of  measured  and  simulated 
dose-efficiency  (D/Q)  relations  in  units  of 
[rad(CaF2)/coulomb  @  1  m]  as  functions  of  voltage 
[MV]  and  at  0°,  20°,  and  50°. 


Shot# 

D/Q 

0° 

20° 

50° 

302 

(C) 

Sim 

190  V2  09 

158  V2'23 

15  V2'75 

Exp 

185  V2  00 

145  V2'20 

16  V2'75 

317 

(Ta) 

Sim 

445  V2  84 

339  V2  88 

117  V3'02 

Exp 

690  V2'5 

450  V2'7 

145  V2-9 

IV.  DISCUSSION 

One  potential  source  of  error  is  the  diode  voltage;  using 
Eq.  (1)  for  voltage  and  an  V2  scaling  for  D/Q,  a  5%  error 
in  measured  current  results  in  a  15%  error  in  the  predicted 
dose  from  Eq.  (3).  This  emphasizes  the  importance  of  an 
independent  measure  of  the  voltage.  Other  sources  of 
error  include  TLD  calibration,  room  scatter,  non- 
equilibrium  electron  incidence  angles  on  the  anode,  and 
ion  current  in  the  diode;  these  will  be  considered  in  future 
work. 
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